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ABSTRACT 


The  behavior  of  the  a-p  tracker  or  filter  which  is  found  in  track- 
while-sc&n  radars  was  investigated  under  various  situatior^.  First 
the  target’s  motion  was  translated  inu>  signals  which  excite  the  filter. 
It  was  found  that  these  signals  could  be  characterized  in  terms  of  tar¬ 
get  performance  capabilities  such  as  velocity,  acceleration,  and 
structural  loading.  Neglecting  false  aianns  and  fades,  the  probability 
of  breaking  urac''  was  found  to  depend  on  the  scan  time,  the  measure¬ 
ment  noise,  the  target  trajectory,  and  the  parameters  a  and  jS.  The 
tracking  errors  could  be  categorized  into  variances  ard  mean  errors 
which  were  investigated  more  fully.  Using  a  postulated  strategy,  the 
effects  of  false  alarms  and  fades  were  found  on  die  variances  and 
mean  errors. 
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BEHAVIOR  OF  a-fi  TRACKER  FOR  MANEUVERING  TARGETS  UNDER  NOISE, 
FALSE  TARGET,  AND  FADE  CONDITIONS 

1.0  INTRODUCTION 

The  subjects  of  automatic  detection  and  tracking  have  received  a  considerable  amount 
of  attention  m  the  last  several  years.  For  instance,  the  NTDS,  MTDS,  and  SPS-33  systems 
are  operational,  and  the  Gilfillan  Company  and  Applied  Physics  Laboratory  (APL)  systems 
have  been  {Hoposed  for  the  SPS-48.  These  systems  have  varying  degrees  of  automation, 
and  with  the  exception  of  NTDS,  have  the  capability  of  automatic  detection  and  tracking 
of  targets.  The  NTbd  sy.stp.n.  vthi'e  having  tracking  equations  and  correlation  regions,  de¬ 
pends  on  a  manual  updating  of  tracks. 

A  system  which  performs  autojnatit  tracking  (not  detection)  is  presently  being  in¬ 
vestigated.  The  automatic  tracking  system  (.ATS)  operates  as  follows.  Tracks  are  initiated 
either  by  an  operator  or  possibly  automatically  if  clutter  is  not  present.  After  initializa¬ 
tion,  an  a-8  tracker  makes  smoothed  esimates  of  a  target’s  position  and  velocity  and  uses 
these  estimates  to  predict  the  target  position  at  the  next  update  time.  A  correlation  region 
is  c-entered  at  tlie  predicted  position  and  an  automatic  detection  is  made  in  the  correlation 
region.  The  advantage  of  this  system  (in  compariscn  to  other  systems)  is  that  there  is  a 
high  a  priori  probability  of  the  target  being  in  the  snwU  correlation  region,  and  conse¬ 
quently,  a  much  hi^er  false-alarm  rate  than  usual  L>  tolerable.  In  other  words,  this  pro¬ 
cedure  avoids  system  saturation  due  to  the  large  lumber  of  false  alarms  which  woiild  be 
inherent  in  a  totally  automatic  detection  .system  It  is  believed  that  a  system  of  this  type 
will  be  able  to  track  targets  through  some  degree  of  clutter.. 

The  purpose  of  this  report  is  to  inve.siipat  ?  tie  performance  of  the  tracker  when 
targets  are  maneuvering  and  the  video  processor  i»  plagued  with  such  difficulties  as  meas¬ 
urement  noise,  false  targets,  and  fading  problems.  In  a  mote  general  sense  one  hopes  to 
learn  what  information  is  required  from  tlic  radar  to  track  adequately. 

Section  2.0  investigates  the  relationship  between  the  characteristics  of  maneuvering 
targets  and  the  signals  as  seen  by  the  radar  Section  3.0  studies  the  errors  as.soriatcd  with 
a-0  trackers  caused  by  maneuvering  targets  and  the  measurement  noise.  In  addition  the 
probability  of  breaking  track  i.  discussed.  Sections  4.0  ard  5.0  investigate  the  perform¬ 
ance  limitotions  due  to  false  targets  and  Cades,  resjicctively.  A  discussion  of  the  re.s-ilts  i.s 
given  in  Section  6.0 

ZO  TRANSLATION  OF  TARGET  CKAR.ACTERIS  flCS  INTO  RADAR 

REQUIREMENTS 

In  describing  the  performance  or  design  of  any  system,  it  is  imixirtant  to  know  the 
characteristics  of  the  signals  which  excite  the  sv'stem.  The  disrassion  given  below  i.s  con¬ 
cerned  with  the  problem  of  obtaining  these  characteristics. 
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In  exanuning  target  irajectories  one  finds  that  there  exists  an  uncountable  number  of 
possible  trajectories  and  that  they  are  random  in  character.  .4i though  in  theory  each  tra¬ 
jectory  co'-ld  be  used  individually  to  investigate  the  system  performance,  it  does  become 
difficult,  due  to  ^he  vast  number  of  trajecto>ies,  to  assimilate  and  use  the  large  amount  of 
inforiTiation.  The  two  procedures  that  are  usually  used  to  circumvent  this  problem  are  as 
follows.  The  first  procesl'ire  is  concerned  with  constructing  a  probabilistic  description  of 
th''  trajectories  and  the  second  proced'ure  is  concerned  with  the  construction  of  a  small  set 
of  trajectories  which  are  typical  in  their  frequency  content.,  Tfie  latter  of  these  two  ap- 
f  proaches  is  considered  in  this  section. 

i 

\  The  discussion  begins  by  making  the  following  a  priori  assumption:  All  taigct  trajec¬ 

tories  can  be  approximated  by  either  straight-iint,  constant-altitude  trajectories,  constant-^ 

■  turns,  or  changes  i'’  altitude  a.s50ciated  with  the  ixrev'ious  two  maneuvers  over  at  least  a  seg¬ 

ment  of  tile  trajectory..  The  number  of  trajectories  is  further  limited  due  to  the  physical 
laws  the  target  mast  obey  according  to  its  design.  In  the  case  of  aircraft,  these  limitations 
can  be  characterized  in  terms  of  velocity,  acceleration,  and  structural  loading  g.  As  a  last 
consideration  the  target  performance  must  be  translated  into  signals  that  the  radar  meas- 
.  ures,  which  are  >fange  and  azimuth.  l‘he  limited  set  of  trajectories  described  above  is  now 

analyzed. 


2.1  Straight-Line,  Constant-AUituile  Targets 

The  geometry  of  the  target  trajectory  is  shown  in  Pig.  2.1,  where  .t  and  y  are  the 
rectangular  coordinates,  r  and  $  the  range  and  azimuth,  and  d,  v,  and  a  are  the  heading, 
velocity,  and  acceleration  of  the  target,  respectively. 


The  foJ’-jwing  quantities  arc  defined. 


v/x-  +  y2 

(2.1) 

tan'”  — 

(2.2) 

X 


0 
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=  V  cos  0 

(2.3) 

=  V  sin  0  . 

(2.4) 

Taking  the  derivatives  of  Eqs.  (2.1)  and  (2.2)  ai’-d  using  Eqs.  (2.3)  and  (2.4)  yield  the 
range  and  azimuth  velocities 


I’fl  =  V  cos  (0  -  <&) 


(2.5) 


and 


e  =  -  -  sin  (0  -  0) . 


(2.6) 


Taking  the  second  derivatives  gives  the  range  and  azimuth  accelerations 

Ojt  =  ~  sin2  (0  ■  <6)  +  a  cos  (0  -  0) 

2 

ao  =  -  sm  (0  -  •S)  +  2  ^  sin  (0  -  0)  cos  (0  -  0)  .■ 
r  r* 


(2.7) 

(2.8) 


If  one  computes  the  total  acceleration  a-p  of  the  target,  one  finds 

Op  =  (flR  ~rvo)iii  +  (rc£)  +2vfiVQ)ig  - 

Substituting  Eqs..  (2.5)  tnrou^  {2.S'  into  Eq.,  (2.9)  yields 

ap  =  acos(0-0)/u  +  a  sin  (0-0)10  . 


The  magnitude  of  ap  is 


joji  -  a  . 


(2.9) 


(2.10) 

(2.a) 


It  is  inU'resting  to  note  that  the  range  and  azimuth  accelerations  as  seen  by  the  radar  can 
be  quite  different  from  the  acceleration  of  the  target.  In  fact  the  n»dar  sees  an  accelera¬ 
tion  just  due  to  the  motion  of  the  target..  The  effect  of  this  apparent  acceleration  is  next 
investigated. 


The  maximum  magnitude  of  the  range  or  azimuth  acceleration  Is  proportional  to 


'^max 


and  the  g y  of  the  timn  Is  defined  to  be 

^?.v  = 


-max 

g 


(2.12) 


(2.13) 


where  g  is  the  gravitational  constant. 
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The  effect  of  this  acccSer  .tion  is  shown  in  Fig.  2.2.  This  curve  shows  that  at  the  near 
tinges  these  apparent  accelerations  coufJ  be  troublesome.;  However,  most  isearch  radars  are 
not  concerned  with  such  short  ranges,  and  these  effects  can  usually  be  ignored.  In  sum¬ 
mary,  for  the  far  regions  one  finds  that  the  range  acceleration  depends  only  on  the  target’s 
acceleration  and  on  us  angle  with  respect  to  the  radar,  'fhe  azimuth  acceleration,  in  addi¬ 
tion  to  depending  on  these  quantities,  al.so  depends  upon  range. 


Fig.  2.2— R^nge  acceleration  at  close  ranges  for  straight- 
line,  constant-altitude,  and  velocity  target  trajectories 


2.2  'I\im  Analysis  in  Terms  of  Circular  Motion 

A  target  moving  in  a  circular  path  with  a  constant  velocity  is  used  to  approzimate  a 
turn.  The  geometry  is  shown  in  Fig.  2.3. 


The  law  of  cosiives  for  the  tnangle  is  written 


1-2  =  p2  .  /j2  -  2pR  cos 


and  can  be  rewritten  as 


r 

7^^ 


2pR 

(p2+R2) 


cos^  . 


(2.14) 


(2.15) 


Recall  the  following  expansion. 
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Fig.  2.3— Gctxn«Uy  of  »  ccnsUnt-g  turn 


the  ran^  srceleration  becomes 
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°R 


i;2 

—  COS  2<^  +  0  . 


for  p  «  R  ... 


(2.20) 


In  observing  Eq.  (2,20)  one  finds  the  same  troublesome  apparent  accelerations  as  found  for 
the  constant-altitude,  straight-line  ta^t  trajectories.  Even  though  the  amplitude  of  the 
harmonies  in  range  can  be  rather  small  at  close  ranges,  the  acceleration  of  the  target  as 
seen  by  the  radar  can  be  quite  large  over  at  least  a  short  interval  of  time.  This  effect  is 
similar  to  a  harmonic  approximation  of  a  square  wave.  Considering  the  first  few  harmonics 
which  are  all  reasonably  small  compared  to  the  fundamental  one,  we  see  tliat  rapid  changes 
in  slope  can  occur  over  small  intervals  of  time. 

The  analysis  of  the  azimuth  measurement  for  a  target  moving  in  a  circle  closely  paral¬ 
lels  the  ran^  measurement  analysis.  One  begins  by  writing  the  law  of  sines  for  the  triangle 
shown  in  Pig.  2.3, 


sin  (0  -9/^)  sin  d 

o  r 


(2.21) 


Solving  for  3,  one  finds 


d  =  Off  +  sin“*  smdj 


(2.22) 


Expansion  into  a  series  yields 


0  I  o 

0  =  Off  +  —  sin  p  -  I  — )  sin^  d  +  0  -. 
"  r  \  r  /  ~ 


(2.23) 


For  the  case  p/r  «  1,  Eq.  (2.18)  car.  be  substituted  into  Eq.  (2.23),  yield-ng 


8.0^.  *  0. 
*'  R  -  p  cos 


(2.24) 


Using  a  dhision  process.  Eq.  (2.24)  can  be  rewritten  as 


0  =  0 

The  azimuth  acceleration  is  computed  to  be 


ff  ^  -  sin  or  +  ^2  sin  2ot  +  0  , 


(2.25) 


1  / \  . 


2/u2\ 


<2.26J 


.Again  one  finds  that  if  p  «  R  and  if  R  is  large,  the  azimuth  var.aticn  :s  nearly  a  sinusoid. 
If  these  conditions  are  not  niet,  the  harmonics  can  cause  a  considerabl  .•  azimuth  deviation 
over  a  small  intcA-al  of  time.. 


It  is  instructive  to  show  an  example  of  tnc  range  suid  azimuth  measurements  and  their 
derivatives  for  a  circular-motion  target  for  two  different  ranges.  The  target  was  moved  in 


NHL  REPORT  7434 


7 


a  circle  at  a  constant  altitude  of  1000  ft,  the  center  of  firvature  was  located  at  an  a;i!nuth 
of  4C'’,  and  ranges  of  14,140  ft  or  141,400  ft  were  used.  The  target  velocity  was  1200 
ft/sec  and  the  gs  of  the  turn  was  6.  The  range  and  azimuth  variations  with  their  derfva- 
tivos  are  shown  in  Figs.  2.4  and  2.5.  Observing  these  curves  r>ne  finds  that  at  the  large 
ranges  the  varations  are  nearly  sinusoidal.  At  the  shorter  distances  it  is  found  that  the 
signals  are  rich  in  harmonics  and  large  acc*e!erations  can  result.  In  addition  it  is  found  that 
the  azimuth  variations  are  quite  small  for  the  far  ranges  in  comparison  to  the  near  ones. 


2.3  Accelerating,  Ciicuur-MotkHt  Targets 

It  is  first  assumed  that  the  target  is  operiing  in  the  far  region  such  that  the  range  and 
azimuUi  variations  can  be  assumed  to  be  sinusoids.  The  target  is  assumed  to  be  moving  in 
a  circular  path  except  that  now  it  may  have  a  tangential  component  of  acceleration.  This 
tangential  acceleration  is  expressed  in  terms  of  angular  acce'eration.  and  is  assumed  to  be 
sinusoidal;  i.e.. 


00  =  —  sill  .  (2.27) 

P 

Integrating  twice  to  find  the  phase,  it  Is  then  placed  into  Eq.  (2.18),  yielding 

sin  >  (2.28) 

p^m  ' 

tn 


V 

wo  -  6  =  —  . 

P 

.\  well-known  Foaner  series  expansion  exists  for  this  waveform  in  terms  of  Elessel  functions 
(1); 


r  =  R  +  p  CCS  (wof 


s 


where 


f  Q  \  /  a  \ 

r  =  R  *  pjn  j - O  1  cos  Unt  +  2pJo  i - »  1  cos  2(CJl)  -  COm  )t 

-  2pJ2  (  ~'7\cos  2(000 +  +  0.  (2.29) 

As  the  tangential  acceleration  goes  to  zero.  Eq.  (2.29)  reduces  to  that  of  the  constant- 
velocity,  circular-motion  target 


r  =  i?  +  pcoscoQt-  (2.30) 

The  tangential  acceleration  in  a  turn  has  the  effect  of  spreading  the  frequency  !.,!ec- 
tra.  The  above  analyss  is  only  valid  for  small  values  of  (a/p^^)  because  of  physical  con- 
sidoetions. 

It  is  probably  while  commenting  cn  targets  changing  altitude..  If  one  allows  no 
azimuth  changes,  the  anal3rsis  so  far  applies  to  th.s  problem  by  merely  exchanging  the  word 
elevation  for  azimuth.  .\n  analysis  of  target  characteristics  considering  a'l  three  measure¬ 
ments  is  not  attempted. 
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Fig.  2. -J— Comparison  of  range  and  derivatives  for  circular-motion 
targets  at  two  langes 


2.4  Frequency  Re^once  Characteristics 

The  discussion  is  begun  by  considering  only  consttiiit-veivv.«ty,  circuiar-motion  targets. 
In  these  cases  only  the  fundamental  sinusoid  is  present.  Thus 


r  = 

R  -  p  cos  pt 

(2.31) 

0  = 

6ft  +  —  sin  0* 
h 

(2.32) 

0  = 

V 

P 

(2.33) 

t:2 

(2.31) 

0  = 

sis'  ■ 

pg«JwgMbn;^M*3Sjmaag?^g»^** 
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Fig.  2.5--Conip«nsor  cf  srimuth  and  derivatives  for  circular-motion 
ur^ts  at  two  ranges 


I 


i 

? 


Subst’futjng  Eq.  12.34}  into  (2.31),  (2.32),  and  (2.33),  one  obtains 


(2.35) 

(2.36) 


Obs€3ving  Eqs.  (2.35)  and  (2.36)  one  fintLs  that  ♦he  amplitude  frequency  characteristics  of 
range  .^-epend  only  on  the  aircraft  parameters  of  the  .structural  loading  g  and  velocity..  The 
azimuth  in  addition  depends  on  range  The  amplitude  Is  plotted  vs  frequency  for  various 
velocities  and  for  sharpness  of  the  tom  indicated  by  in  Fig.  2  6.  One  finds  that  as  gy 
for  fixed  v  is  increased,  the  amplitude  decreases  and  the  frequency  increases.  .\s  the  ve¬ 
locity  f  '>T  fixed  g\  is  increased,  the  amplitude  increases  and  frequency  decreases. 
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.-U!  example  of  the  trajectory  of  a  high-performance  aircraft  performing  “dogfight” 
type  mareuvers  is  given  ir<  Fig.  2.7.,  The  range  and  azimuth  variations  with  time  are  given 
in  Figs.  2.S  and  2  9  for  this  maneuver.  This  e.xampie  again  shows  that  the  frequencies  en- 
cc-untered  in  even  high-perlormance  targets  are  quite  low. 


2.5  Discussion  of  Results 

The  purpose  of  this  section  is  to  translate  the  target  characteristics  into  Jgnals  on 
which  the  radar  must  track.  It  was  found  »hat  the  three  target  parameters  g,  v,  and  a  can 
be  used  to  characterize  the  limitations  on  the  received  signals  in  terms  of  amplitude  and 
frequency  characteristics.  If  the  near  effects  could  be  ignored,  the  s^als  set  up  by  the 
target  were  foimd  to  be  quite  low  in  freqrrency;  typical  values  may  be  in  the  ordt.-  of  0.2 
Hz  or  lowCT.  Howeva",  if  the  target  has  a  tangential  component  of  acceleration  as  well  as 
a  turning  motion,  hi^er  harmonics  in  the  radar  signal  may  be  observed  in  which  the  hip¬ 
est  sgnificant  observjble  frequency  may  be  0.6  to  0.8  Hz. 

It  is  assumed  threu^out  the  remainder  of  this  discussion  that  only  the  polar  coor¬ 
dinate  system  will  be  used  for  tracking  and  that  the  near  effects  for  all  practical  purposes 
can  be  ignored. 
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Fig.  2.8— Tinic  history  of  rangt  for  a  targvt  performing  iramni- 
rers  in  a  iim't«i  spacn.  rargt  about  .I.S  naut.  mi. 
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T!Mt 

Fig.  2.9— Tine  history  of  azimuth  for  a  target  maneuvering 
ui  a  limited  space;  range  about  35  naut.  mi. 


3.0  OPERATION  OF  TRACKING  EQUATIONS  EXCLUDING  FALSE  TARGET 
AND  F.ADING  CONDITIONS 

In  grnoral  the  tracking  equations  niusl  perform  three  functions.  They  must  compute 
the  velocity,  .smooth  both  Uie  measured  posiuon  and  computed  velocity,  and  fin?Jly  pre-. 
diet  a  new  coordinate  on  which  to  center  the  correiation  region.  The  a-0  tracking  equa¬ 
tions  are  first  defined  and  the  errors  in  both  the  mean  and  v'ariaiice  of  the  predicted  coor¬ 
dinates  are  irvc  -»igated.  Finally  the  probability  of  breaking  track  is  discussed. 


3.1  Definition  of  Tracking  Equations 

The  well-known  a-p  tracking  equations  are  defined  as 


y.V  _  v'-l 

+  vf-’r 

(3.1) 

4 

< 

II 

a(Af  ) 

(3.2) 

if  =  if-’ 

7‘ 

(3.3) 

where 

Xp  =  predicted  position 
Xj  =  smoothed  position 
V,  =  smoothed  velocity 
X„  =  measured  position 
a.p  =  filter  na’ameters 
T  =  sampling  time 
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For  purposes  of  analysis,  these  equations  are  mo>Te  usable  if  one  subsniutes  Eq.;  {3.1)  into 
Eqs.  (3.2)  and  (3.3).  yielding 


and 


[(1  -a) 
! 

!  -S/T 


{i-^?j[vij  i^/rj 


,  S*l 

4 

Pi 


=  [i 


'  v  i 


.3.4) 


(3  5) 


For  convenience,  Eqs.  (3.4)  and  (3  5)  are  put  in  block  diagram  form  in  Fig.  3.1  and  the 
following  matrixes  and  vectors  are  defined: 

HI -a)  (l-a)rl 

A=  I  i 

L-jJ/T  (l-i?)  J 


B 

C 


! "  1 

[iJ/rj 

[l  T-j 


X 


.At  this  point  the  system  has  been  reduced  into  a  more  tractable  form  for  analysis. 


Fig.  3.1— Block  dtagram  of  an  o-p  irackcr 


3.2  Stability  .Analysis 

Since  the  system  is  in  the  form  of  differenc?  equations,  it  is  proposed  to  investigate 
the  stability  through  ti’e  use  of  the  Z-Uansform  (2».  Taking  the  Z-transform  of  Eq.  (3.4) 
results  in 


X(Z)  =  .\Z-i  .XfZl  +  BA„,(Zl. 


(3.6) 
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1  ne  system  transfer  function  is  doftned  as 


G(Z) 


X(Z) 


(/ -.4Z-M"’S  . 


For  tiie  sj>ecific  rv-stem  in  question,  Eq.  (3.7)  becomes 


G.,(Z‘  = 

A's(Z) 

oZ(z  -  ] 

\  g  / 

-V,„(Z) 

Z2  Det  (/-  .4Z-M'' 

G,{Z)  = 

VsxZ) 

i  ziz-i) 

-V„.(Zt 

~/^6ei[T~AZ^\^- 

where 

Z~  Del  (/-.4Z-’]  =  Z"  -  Z(2-q  -p>  +  <1  -a),. 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


The  stability  of  the  system  depends  on  the  case  for  which  at  least  one  root  of  the 
Z^  Del  (/  -  .4Z-1 1  lies  outside  the  u.nit  circle.  Several  typical  root  loci  are  shown  in  Figs. 
3.2  and  3.3,  .-Vs  can  be  seen,  the  parameter  a  seems  to  control  the  radius  of  the  loci  end 
P  controls  the  amount  the  pole  has  moved.  For  cases  in  which  a  >  1,  the  loci  never  eaves 
the  real  a.\is.  The  stability  cntenon  can  be  found  by  setting  the  most  negative  root  equal 
to  -  1, 


Z  =  -! 


(2  -g  -^) 
0 


~  x'Ta  4  ^)-  -  A3  . 


Solving,  one  finds  that 


(3.11) 


3  %  A  -  2a  ,  (.3.12) 

In  addiiion  to  the  stability  criterion,  ?.  retc  in  Eq.  (3.7)  can  lie  outside  the  unit  circle. 
This  situation  will  not  Oe  allowed  since  u  usually  leads  to  undesirable  system  rt^ponse  (2'. 
The  condition  one  must  satisiy  is 


P  <  2u  .  (3.13) 

The  permissible  values  of  a  and  3  are  shown  :n  Fig.  3.4.  Further  restrictions  will  be-  placed 
on  this  region  later. 


3.3  Frequency  Respot<se  Characteristics 

The  frequency  respon-^  of  the  filter  can  be  found  by  placing  Z  =  into  Eqs.  (3.8) 
and  1 3.9).  The  magnitude  and  pha.se  characteristics  of  O^iZ)  and  G,  (Z)  are  shown  in  Fi^. 
3.5,  3.6,  and  3.7. 


Observing  these  figures  Oi.e  finds  that  a  should  never  be  larger  than  one.  For  q  <  1. 
g  seems  to  control  the  bandwidth  of  the  low-pass  filler  and  3  has  more  control  over  the 
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Fi?.  -3.S--Amp:ilude  and  phase  rf-.aractenstits  of  tJie  tracking  equations 

fora  =  0.9 


damping.  In  fact  ^  should  be  some  what  smaller  than  rs  such  that  resonant  spikes  do  .not 
occur,  Obserting  Rg.  3.5.  one  finds  that  the  measured  position  is  nearly  unsmoothed  and 
that  Gi-fZ)  Is  in  the  form  of  a  differentiator.  In  general  a  and  P  ^tould  be  set  such  that 
Ox(Z)  passes  the  highest  expected  freq»jency  andistorted,  and  G^iZ)  acts  as  a  good  dif- 
ferentiatoc  up  to  the  highest  expected  frequency. 


3.4  Noise  Characteiistics  of  Filter 

The  discussion  is  begun  l:<y  '  aking  the  foliowing  chanee  of  variable; 

iX  kj/.  {3.14} 

Equations  (3.4)  and  (3.5)  describing  the  filter  can  be  rewritten  as 


rx.i-'- 

"(1  -  a) 

(1- 

■a.i] 

V  '  -1  /V  ~  V  T  .V 

X*  !  i  O  -iXml 

i 

1 

1 

i  1 

'•  !  i 

(3.15) 

iLXj 

L 

(I  - 

■m\ 

LXj  i_3j 

[1  1] 

(C.16) 

ij 
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Fig.  3.0— Amplitude  and  phase  chan^HeniUcs  of  tracking  equations  for 

o=  O.I 


where  now 


[(l-a)  (1-0)1 

j  I 

[  -a  (l-^)J 


Tal 
B’  =  l  i 

L^J 


-=[1  1] 


The  mean  and  rovariance  equations  are  defin.'d  as  '3) 


[X]^  -  VX'-*  +  B'X^ 

r't 

p.v  =  a’P'-'A’T'  +  BV-S' 


where  [a']-'  =  £'[x'']  ,  the  expected  vaJuo. 
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a<T/4  iiiT/2  u»T/4  wT/Z 


wT/4  wT/2  Ci;T/4  «T/2 

Fi^,  3.7~Ainpl:t  jde  and  phase  cha.'ac*»TisUcs  of  Uie  tracking  equations  for 

a  •=  1.3 


r^ll  ^12 1 

|_^2!  P22J 

p-,l  =  E[{Xs-XsHXs-Xs)] 

P\2  =  £[(A', -.Yj(A.V-Zi.V  J  =  P21 

P-ll  =  £[(AA-- AX1(AA'-aJ)j 

o.^  =  siandard  deviation  of  measurement  error  . 

The  mean  equations  are  just  the  response  of  the  system  to  a  deterministic  input  which  ’Ai}! 
be  studied  later.  The  covariance  equation  becomes 
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Tlie  variance  of  the  predicted  positit  n  is  easily  computed  in  terms  of  the  variances  ^ven 
m  Eq.  (3.19): 

=  I’ll  +  2Pi2  *  ^22  >:  (3.20) 


Ife.noring  the  transients  during  the  time  the  track  is  being  established,  the  steady  state 
solution  is  found  by  recursively  solving  Eq.  (3.19).  The  predicted  noise  power  is  nonrial- 
ized. 


J3.21) 


and  plotted  as  a  function  of  or  and  0  in  Fig.  3.8.  Observing  Fig.  3.8  one  finds  tliat  the  in¬ 
put  noise  power  can  be  reduced  or  filtered  if  a  and  p  are  appropriately  adjusted.  One  also 
finds  that  the  parameter  P  greatly  iffects  the  predicted  noise  power.  This  effect  should  be 
expected  since  differentiated  noise  is  quite  ncisy  and  p  affects  this  quantity.-  In  fact  the 
variance  of  the  predicted  position  can  be  larger  than  the  mput  noise  because  o  f  differen¬ 
tiation. 


ALPHA,  a 


Fig.  3.8- Predict^  posiUon  noise  power  as  a  function  of 
aj 


An  interesting  phenemenon  occurs  for  sm:ill  values  o*"  a  and  large  taiues  of  3.  The 
noise  p>ower  increases  sharply.  The  reason  for  this  may  most  easily  be  seen  by  look'ng  at 
the  frequency  responses  in  Fig.  3.6.  The  large  resonant  peaks  in  the  resfionse  aJlow  a  lot 
of  noise  to  come  throu^.  In  addition  a  non-minimal  phase  conditio;,  can  be  reached  and 
the  noise  then  increases  rapidly.. 

The  noise  power  is  not  a  function  of  the  sampling  time.  .As  the  sampling  time  de¬ 
creases  the  variance  of  the  velocity  increases;  however,  this  effect  is  canceled  out  by  having 
to  predict  over  a  shorter  interval  of  time. 
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3.5  OeUrministic  Eiror  .Analysis 

Various  sinusoidal  waveforms  are  passed  throu^  the  filter  and  the  root  mean  square 
error  between  the  predicted  and  true  positions  is 


where 

X*p  =  sin  +  +  IGfOto)!  sin  (/ajT +  dGj,) . 

given  that 

=  sin  ((i  +  DcoT)  . 

(the  ampJiU'de  of  the  sinusoid  equals  one  for  normalization) 

"Hie  rms  error  is  plotted  vs  the  ratio  of  the  frequency  f  of  the  sinusoid  f  to  the  sampling  fre¬ 
quency  ff  for  various  values  of  a  and  j?  in  Figs.,  3.9  throu^  3.11. 

Jpon  inspecting  these  curves  one  finds  that  if  eitha-  the  waveform  is  sampled  faster 
or  thi  frequency  of  the  waveform  is  lessened,  then  the  error  decreases.  In  addition  the 
lightlv  filtered  systems  have  fewer  errors  in  the  means  than  the  heavily  filtered  systems. 
The  regions  of  the  curves  where  the  error  is  near  unity  only  indicate  that  one  is  not 
sampling  at  a  sufficient  rate. 


001  03  IOC 

FBEOUENCY  RATIO,  f/ty 

Fij.  3.9— IVt^nainHlic  prediction  error  a«  a  functkm 
of  f/f-p,  a,  p 
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Fig.  3.10— Deterministic  prediction  error  as  a  func¬ 
tion  of  f/fT,  a,  ^ 


3.f  Ptobability  of  Breaking  Track 

A  correlati(Mi  region,  as  shown  in  Fig.  3.12,  is  defined  to  be  a  region  of  space  cen¬ 
tered  about  the  predicted  value.  If  iA’m  ~  A’p  I  <  y,  then  the  measured  value  is  said  tc<  lie 
in  the  ooneli.tion  r^on. 

The  random  variables  and  Xp  are  gaussian-distributed  random  variables  as  snown 
in  Fig.  3.13.  The  random  variable  to  is  deHned  as 

u  =  Xn  -  Xp.  (3.23) 

Since  A'p  and  Xm  are  gaussian,  cj  is  also  gausaan.  The  mean  and  variance  of  u  are 


Cii  —  Xf^  Xp 

(3.24) 

9  O  O 

of-  =  Ot  . 

u-  -vp  Am 

',3.25) 

The  probability  cf  breaking  track*  at  a  particular  sample  instance  is  defined  as  the 
[Hobability  that  |X„  -  .Ypl  >  y. 


*Thit  tmc.mes  that  once  the  target  is  outside  the  correlation  regr-m  it  will  remain  outside  oi  it  for  aO  time. 
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Fig.  3  ■  I— Deterministic  prediction  error  as  a  function 
of  f/ff.  a,  C- 


CCRPELATION 

REGION 


1  X 

-r  Xp  Xp+r 


Fig.  3. 1 2— Corre^ition  region 


Xj-r  x„  Xt,  x„  x,+r 


Fig.  ^.IS-Profaat-ility  densitier  of  Xp  and  define  i 
over  correlation  region 


CO  Xfn  Xp  —  f{<X,^\T,Xni) 

oL  =  4p  =  Aa..i^Oxm)  ■ 

Observing  Eq.  (3.8)  one  mi^t  infer  that  the  real  problem  in  minimizing  the  protebility  of 
breaking  trac.  for  a  given  T  and  measurement  noise  is  that  the  target  trajectory  is  so  un¬ 
predictable.  Digressing  for  a  moment,  we  discuss  below  several  possible  approaches  to  filter 
design. 


Some  filters  assume  that  the  tablet  trajectorj'  is  known  and  adjust  q  a^nd  3  in  an  adap¬ 
tive  way  such  as  to  minimize  the  error,  i.?.,  the  Kalman  filter  (4).,  Some  adaptive  filters 
try  to  estimate  the  future  bandwidth  requirement  (how  fast  X„,  is  varying!  and  adjust  a 
and  8  to  minimize  the  error.  OthCT  filters  assume  a  given  trajectory  and  mininiirte  a  cost 
function  such  as  to  determine  8  in  terms  of  q  (5).  Other  filters  use  fixed  c  and  p,  ad¬ 
justed  such  that  the  error  is  minimized  for  the  worst -case  target  trajectory;  the  error  ’S 
bounded  in  this  manner.  In  this  author opinion  mmy  of  the  types  of  filters  discussed 
are  not  well  suited  to  track-while-scan  radar  operation.  Most  of  the  difficulty  se‘=ms  to 
arise  from  assuming  a  priori  target  trajectories  or  making  fat-ire  predictions  based  upon  loo 
little  'pformafion.  .-X  filter  design  will  be  discussed  in  a  sub-iequenl  report. 


3.7  Some  Examples 

The  a-2  tracker  is  used  to  track  the  hi^i>  maneuverable  targets  described  in  Secuon 
2.4.  The  true  and  predicted  range  and  azimuth  are  plotted  vs  time  m  Figs.  3.14  through 
3.17.  The  correlation  region  was  made  sufficiently  l<uge  .such  that  the  track  would  not  oe 
broken.  The  sampling  time  is  4  sec.  Two  different  valuer,  of  a  and  3  were  used.  The 
standard  deviation  of  range  and  azimuth  ncise  are  250  ft  and  0.52  degrees,  respectively. 
The  fiitered-noise  standard  deviation  a^p  is  shown  with  ihc  predicted  position.  The  mcas- 
uiement  noise  of  the  true  trajectory  is  not  shown  since  it  tends  to  make  Uie  graph  difficult 
to  read. 

Observing  these  figures  we  find  that  the  lightly  filtered  case  had  ic-s.s  error  in  the  mean 
but  a  larger  staiidard  deviation  Oyp  than  the  more  heavily  filtered  case 

The  same  tracking  situation  as  found  in  Fig.  3.14  is  again  shown  in  Fig.  3.18  except 
that  the  sampling  time  is  reduced  to  1  sec.  One  finds  that  the  mean  of  the  predicted 
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Fig.  3. 1  4— CompArison  of  true  and  predicted  signals  for  a  highly  maneuverable  target 


•64)00- 


Fig.  3. 1 5— Comparison  of  true  and  predicted  sisals  for  a  hig£.iy  maneuverable  target 
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TiKe  ;**c) 


Fig.  3. Ifi— Comparison  of  Srie  aid  priHiicted  azimuths  for 
a  hi^ly  maneuver;  ble  target 


position  tracks  very  close  to  the  true  trajectory  (the  standard  deviation  of  the  noise  is  the 
same  as  found  before).  This  agroes  with  the  theor\’  in  that  a  decrease  in  samj  ling  time  re¬ 
sults  in  improved  system  performance.  The  same  effect  could  be  achieved  if  t!:e  target 
did  not  maneuver  quite  so  violently  such  that  the  frequency  content  of  the  signal  was  less. 
This  example  shows  how  critical  it  is  to  know  the  1  requency  content  of  the  target  trajec¬ 
tory  in  order  to  design  a  good  tracking  system. 


3.8  rhscus-sion  of  Results 

The  behavior  of  the  o-J  tracker  excluding  the  effects  of  false  targets  and  fades  was 
discussed.  Althou^  the  analysis  was  quite  straightforward,  no  reasonable  way  of  adjust¬ 
ing  a-P  to  minimize  the  probability  of  breaking  track  for  a  given  7  was  found.  This  diffi¬ 
culty  arose  frcm  the  fact  that  the  probability  of  breaking  track  dejiended  on  the  target 
trajectorj  and  therefore  the  minima  would  be  different  for  each  trajectory.  .-Xdaptive  and 
worst-cate  design  solutions  were  mentioned.  Several  examples  were  given  to  illustrate  the 
theory. 
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4.0  SYSTE.M  LIMITATION  JUE  TO  CLLITER 

.4s  5ton  in  tht  section  one  could  make  the  probability  of  breaking  track  as 

small  as  one  wishes  by  making  the  correlation  region  larger.  In  reality  this  is  not  the  case 
because  false  targets  cau.sed  by  clutter  or  other  targets  place  a  basic  limitation  on  this  quan¬ 
tity.  This  effect  is  investigated  next.: 


4.1  Average-Target  Strategy 

It  is  not  known,  at  least  to  this  author,  what  constitutes  an  optimum  strategj’  for 
processing  multiple  targets  within  the  correlation  region  nor  how  one  weald  construct  such 
a  procedure.  As  an  alternative  to  searching  for  an  optimum  processor,  a  strategy  is  i>os- 
tulated  and  analyzed.. 

It  is  assumed  that  any  target  present  in  the  com-lation  region  is  an  equally  likely  can- 
didai-*  for  being  the  trtie  target.  In  other  words  no  a  i>riori  knowledge  is  av'ailable  and  the 
measurements  are  indistinguishable..  Since  it  has  been  assumed  that  all  targets  are  equally 
likely,  the  information  from  all  targets  will  be  used..  The  strateg>-  postulated  to  achieve 
this  goal  is  to  average  the  position  of  aU  the  targets  in  the  correlation  region;  i.c.. 


'  (A'  +  ll  ’ 

where 

K  -  number  o;  false  targets  in  correlation  region 
n,  =  zero  mean  uniform  noise  shown  in  Fig.  4.1., 

becomes  the  filter  input  rather  than  as  previously  discussed.  Tne  filter  equa¬ 
tions  can  be  written  as 


(Al'  =  +  B\Xi]^ 

1 

^ — Ihr — '■ 


•'''  =  r-  f  v'l  -V-l 


lApi'  =  cm 


Equations  (4.2)  and  (4.3)  are  comoined.  yielding 

-yT'’  [(\ _ SL^\  _  _SL_\^rvT'“^  [  iL.v  ^  N”* 

\V  A  +  l)  K*ir 

i  .  -I  |L 

!v  !  A  -  V  i  I-  "'-i 

:  *J  !  {K*1)T  \  K*lj  *  i  r(A  +  1)  1 

■■  L  J  i _ I  L  _J 


j  (4.5) 
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i 

p(n.) 

_ _ _ _ -  _ 

- ^ 

Xp-y  Xp  Xp+y 

Fi*  -(.1  -Probability  density  function  for  clutter 


Define  a'  and  ^  as 


a 


(Jf  +  1) 

/r  -  -J-_ 

(ff  +  i)  “ 

By  using  these  definitions,  Eq.,  (4.o)  becomes 


r  ^ 

rA,  1'^'  fd-o')  (l-a')r 

[x,  !  a' 

i  1  ♦  > 

1 

L^J 

L  '=1  J 

(4.6) 

(4.7) 


(4.8) 


!n  this  etjuation.  the  number  of  false  tai^ets  K  controls  the  amount  of  smoothing  and  the 
amount  of  measurement  uncertainty. 


The  solution  for  the  noise  power  of  Eq.  (4.5)  is  not  nearly  as  easy  as  before.;  The 
reason  is  that  It  is  a  random  variable.  Equation  (4.8)  involves  the  product  of  two  random 
variablts  and  in  general  the  probability  density  changes  form  from  one  iteration  to  the 
next.  For  this  reason,  only  special  solutions  will  1  e  found. 


4.2  Noise  Characteristics  for  Given  Clutter  Conditions 

S{)e<ia!  solution-^  to  Fxj,  (4.8)  can  be  obtained  by  assuming  that  K  is  a  deterministic 
sequence.  For  example. 


1  0  0  0  0  1  0  0  0  0  1  0  U, , , 


4  4  0  0  !  4  U  0  4  4  0. , 


The  equation  becomes  a  stochastic  difference  equation  with  time-vuiysf-f  coefficients.  In 
addition,  the  uniform  probability  dciisity  of  the  clutter  is  approximated  as  gaussian  with 
the  variance  set  equal  to  the  variance  of  a  uniform  variable.  The  covariance  equations  are 
found  in  the  same  manner  as  before,  i.e,. 
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.V 

“(l-a')2 

2(l-a')2 

(l-a')2 

^11 

P\2 

= 

-^(1-a) 

(l-Z^Ml-a) 

(l-?'Ml-a) 

P\2 

J*22_ 

^2 

-2^(1-^) 

(1-^)2 

P22 

+ 


+ 


_  „2 

(if  +  1)2  % 


^2 


where 


(4.9] 


a'  =  aHK  +  l) 

^  =mK  +  i) 

=  7/\/3 

=  a  deterministic  sequence  taking  on  thi;  caiue  of  the  number  of  fa^  alarms  at 
^ch  iteration. 


The  foUowin » terms  are  defined: 


oo  = 


«.Y„ 


Ot  = 


(4.10; 


(4.ii; 


The  predicted  noise  piwer  is  plotted  vs  time,  for  various  conditions,  in  Fi^.  4.2  through 
4.9. 


Fig.  4.2— Noise  power  as  a  function  of  time  for  requence 
10000010. . 


samples 


Fi^  4.3— Noise  power  as  a  function  of  time  for  seque?Mx 
40000040. . . 


Fig.  4.4— Noise  power  as  a  function  of  lime  for  sequem-e 
lOOOOOiO..  . 


Fig.  4.5— Noise  power  as  a  function  of  lime  for  srqu-.ice 
40000040 . . . 
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Fig.  4.9— Noise  power  ss  a  function  of  li.Tic  few 
sequence  4-HGOCOOC4  440 . . . 


These  figures  show  that  when  false  targets  appear  periodically  then  a  steady  state  solu¬ 
tion  exists  for  the  noise  power  of  the  predicted  position.  The  noise  power  increases  dur¬ 
ing  the  time  the  clutter  is  present  (it  will  be  shown  later  that  the  no^se  power  saturates  at 
a  given  value  in  extended  clutter!,  and  drops  back  to  its  steady  state  value  after  the  false 
targets  are  removed.  In  all  eases  it  was  found  that  when  the  correlation  n^on  size  became 
larger  as  indicated  by  o^,  the  predicted  noise  power  became  larger  during  the  period  of  the 
false  alarms.  TSie  heavy  smoothing  case  (a  =  0.2.  p  =  0.05)  was  generally  less  noisy  than 
the  light  smoothing  (a  =  0.9,  3  =  0.7)  case.  The  sequences  with  extensive  clutter  (more 
than  one  false  alarm  in  a  row)  were  generally  much  more  noisy  than  the  single  false  alarm 
case. 


If  the  £e«iuences  are  such  that  K  is  the  same  value  for  all  sampling  instances,  the  pre¬ 
dicted  noise  power  saturates  at  a  given  value.  These  results  are  shown  in  Tigs.  4.10  and 
4.11.  fliess  figures  indicate  that  larger  correlation  regions,  more  false  alarms,  and  lifter 
damping  all  result  in  larger  predicted  noise  power. 


4.3  Deterministic  System  Behavior  Under  Clutter  Conditions 

If  one  either  finds  the  mean  value  or  eliminates  the  noise  in  the  system  described  by 
Eq.  (4.8)  one  has  the  following  equation: 


1(1 -o')  (l-a')Tj 

“  - 

•v-i  fo'' 

1  * 

i 

1  (i-3')J 

V 

Iff/T 

(4.12) 


The  system  is  excited  with  a  sinusoid  and  ihe  root  mean  square  error  is  found  in  the  same 
manner  as  described  in  Section  3.5.  The  only  difference  .n  the  solutions  obtained  here  is 
that  a'  and  /J'  vary  from  sample  to  sample  according  to  a  deterministic  sequence  in  fC;  i.e.. 
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where 


* 

a 


tt 

K  +  i 


H’ 


K  =  1100011000110 


(4.13] 


The  e.Tors  are  shown  as  a  function  of  various  false  alarm  sequences  and  parameters  a  and 
.3  in  Figs.  4.12  and  4.13. 


F'.g.  4.12— Mean  value  frequency  response  under 
false  aUrn)  conditions 


Fig.  4.13— Mean  value  frequency  response  under 
false  alarm  conditions 
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It  is  found  that  the  error  in  the  mean  grows  as  the  number  of  false  alarms  increases. 
Of  course,  it  is  again  Tound  that  as  tiie  sampling  frequency  increases  or  the  excitation  fre¬ 
quency  decreases  the  errors  are  less.  Also,  the  more  ii^tly  filtered  system  has  less  error 
than  the  more  heavily  filtered  astern. 


4.4  Probability  of  Breajdng  Track 


In  Section  3.6  the  probability  of  txeaking  track  was  foimd  to  be 


P[break  track] 

In  the  presence  of  clutter  both  iuj  and  0(^ 
track  worse. 


1  1  _ 

2  -  2 

increase,  making  the  probaoility  of  breaking 


It  is  not  always  possible  to  make  the  {Hobability  of  breaking  track  smaller  by  making 
correlation  region  larger  as  can  be  seen  from  the  following  argument.  As  one  increases  the 
size  oi  the  corrclaiion  region,  the  chances  of  receiving  one  or  more  false  alarms  increases. 
.\f  some  point  the  probability  of  Ineaking  track  would  start  increasing  with  increasing  7 
due  to  th**  iai^  enters  set  up  by  the  false  alarms.  'iTius  for  a  given  clutter  environm'  nt, 
the  probabihty  of  breaking  track  can  not  be  made  arbitrarily  small. 


4.3  An  Example 

The  same  tracking  situation  as  shown  in  Fig.  3.16  is  used  except  that  in  this  case  four 
false  alarms  per  correlation  region  are  assumed  to  occur  at  t  =  16,  20,  and  24  sec.  The 
correlation  region  is  assumed  to  be  16  times  larger  than  the  measured  standard  deviation.. 
The  resulting  track  is  shown  in  Fig.  4  14. 

One  finds  that  whei*  clutter  is  present,  the  filter  has  hec/y  smoothing  and  does  not 
negotiate  the  turn  well.  In  other  words  the  mean  of  the  predicted  position  ;s  further  i. 
enor.  In  addition,  the  standard  deviation  oxp  of  the  filtered  noise  grows  diuing  the  time 
the  false  alarms  are  present  due  to  the  added  noise.  Of  course  if  the  heaver  clutter  would 
have  appeared  d  'ring  a  time  that  the  target  was  not  maneuvering  so  violently,  the  errors 
would  hai'e  been  n.udi  less.  This  su^csts  that  the  system  could  be  improved  by  sampliiig 
faster-  There  is,  however,  an  up.aer  limit  to  this  procedure,  and  that  occurs  when  the  clutter 
begins  to  be  correlated. 


4.6  Discussion  of  Results 

An  average-target  strateg>'  was  postulated  and  analyzed.  Because  a  general  solution 
for  the  predicted  noise  power  and  deterministic  tracking  error  is  extremely  difficuit.  only 
special  solutions  were  found.  It  was  fornd  that  the  prediction  noise  increased  when  the 
amount  of  clutter  increased,  when  the  correlation  region  was  made  larger,  and  when  the 
system  used  tighter  filt«ing.  In  addition,  the  deterministic  error  grew  as  the  clutter  levels 
increased,  as  the  system  used  heatier  smoothing,  and  as  the  sampling  frequency  war  dc 
creased.  It  was  argued  that  the  probability  oi  breaking  track  could  not  be  made  arbitrarily 
small  in  the  presence  of  clutter  by  making  the  correlation  region  larger. 


!*»i  njiUBi^TC'lk  uuslHiWflp^^PSiWBHS 
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Fig.  4, 1 4— Comparison  of  inie  and  prvlictcd  >.tni!fi.  for  a 
highly  maneuverabi?  target  in  the  presence  of  ciuite. 


5.0  SYSTE.M  LI.MITATJON  DUE  TO  F.4DES 

Targft  fading  can  be  caused  by  many  factors  of  which  the  multipath  effect  is  possibly 
one  of  the  most  troublescme.  Considering  this  fading  source  only,  one  investigates  the 
parameters  that  effect  the  interval  of  lime  during  which  the  fade  occurs.  The  effects  of 
multipath  fades  on  the  .system  performance  are  analyzed  by  postulating  a  stra'egy  and 
determining  the  characteristics  of  both  the  noise  and  the  dctcrminhstic  errors.,  ft  is  assumed 
that  no  false  targets  are  present.. 

5.1  .Multipath  Fading 

The-  .signal-lo-noise  ratio  at  the  receiver  can  be  computed  from  the  radar  range  enu?- 
tion  oiierating  in  the  noise-limiled  region  (6); 


■  §-  \  =  f 

i  .V  t  "  ( -in  y'KT^Bs  R  -  L 
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where  the  pattern  propagation  facto*s  Ft  and  Fr  are  due  to  the  multipatli  effect.  They 
are  found  by  addiiig  the  electric  field  from  both  the  direct  path  Ea  and  reflected  path  £/> 
as  shown  in  Fig.  5.1; 


Fig.  5.1  Multipath  effect 


These  patterns  have  been  studied  quite  extensively  by  Blake  (7).  Several  antenna  patterns 
reproduced  from  his  work  are  shown  in  Figs..  5.2  through  5.4.  Only  two  of  the  more  im¬ 
portant  parameters,  frequency  and  antenna  hei^t,  which  effect  the  patterns  are  listed  with 
each  figure. 


An  example  is  next  considered.  A  target  is  flown  from  over  the  horizon  toward  the 
radar  at  a  constant  altitude  of  10,000  ft.  The  radar  is  assumed  to  be  a  two-dimensional 
continuous-scanning  search  radar  with  an  antenna  height  of  80  ft  and  a  frequency  of  1300 
MHz.  !f  the  signal-to-noise  ratio  is  above  13  dB  (corres}x>nds  to  3  Pj  =  0.9  and  Pfa  =  10“^> 
a  detection  is  assumed  to  be  made;  below  13  dB  a  fade  is  said  to  occur.  For  the  target 
considered,  the  signal-to-noise  ratio  is  assumed  to  be  17  dB  at  the  renter  of  the  first  lobe 
in  w’hich  the  target  is  detected,  the  lobe  structure  is  computed  using  Blake's  programs.  The 
detection  zones  are  shown  in  Fig.  5.5. 


One  finds  that  both  the  detection  and  fade  zones  -Jeexease  in  size  as  the  range  is  de¬ 
creased.  The  main  -cause  for  this  is  that  the  angle  of  the  lobe  or  angle  between  the  lobes 
multiplied  by  the  range  gives  the  distance  across  the  lobe  or  the  distance  between  the  lobes. 
These  distances  decrease  as  the  range  decreases.  The  fad»*  zone  decreases  more  rapidly  than 
the  detection  zone  with  range,  niis  can  be  attributed  to  an  improved  i  ignal-lo-noise  ratio 
due  to  the  decreased  range.  From  Uiis  simple  example  one  may  estimate  that  fade  zone 
distances  may  range  from  10  naut.  mi.  to  several  thousand  feet. 


The  length  of  time  the  target  fades  is  given  by 
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Fig.  S.4— Aatenna  pattern  fr>r  u  antenna  hei^t  of  20  ft  and  a  frequency  of 

100  MHz 


■TKJ-nM 

..4l 


MO«lZ0 


400XXX) 


SOOa'^ 

RANGE 


soaooo 


DETEC 


-  DETECTION  i 

1 

OETBCTBN 

Dacron 

OETEcflON'j 

1 

> 

fMiM 

fmm 

1 

HOETcCTIOrj 


200,000 


300,000 
FOLDED  range 


Fig.  5.5— Detection  and  fade  zones  vs  ra.ige  for  a  target  at  10,000  ft,  radar  height  80  ft, 
frequency  1300  MHz,  threshold  S/.V  cf  13  dB,  and  S/V  ratio  at  580,000  ft  of  1?  dB 


^  _  Distance  between  ietection  r^ons 
Target’s  radial  velocity 

This  fade  time  is  plotted  ire  radial  velocity  in  Fig.,  5.6  for  a  lO-naut.-mi.;  and  a  10,000-ft 
fade  zone. 


■Jsing  the  same  example,  we  find  that,  at  the  far  .anges  where  the  fade  zones  are 
quite  wide,  the  signal  can  fade  from  about  10  min  at  the  low  vetocitics  to  0.5  min  at  the 
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Fig.  0.6— Signal  tiis  time  vs  radial  velocity 


high  velocities.  At  the  closer  ranges  where  the  fade  zones  are  much  narrower,  the  signal 
can  fade  from  about  s  minute  to  several  seconds  or  lower  at  the  hi'her  velocities. 

It  is  desiraUe  for  later  considrations  to  construct  a  sample  sequence  in  which  a  *‘l” 
denotes  a  fade  and  a  “0”  a  detection.  Ta»gets  are  flown  throu^  the  detection  and  fade 
zones  at  various  velocities.  The  position  of  the  target  is  sampled  every  T  seconds,  accord¬ 
ing  to  tlie  scan  tiine  of  the  radar.  A  1  or  a  0  is  assigned  to  each  sample  instant  according 
to  whe*^her  the  target  is  in  the  fade  or  detection  zone..  Several  sequences  are  given  in  Figs. 
5.7  and  5.8.  The  effects  of  fading  sequences  on  the  noi'"  properties  of  the  filter  are  now 
examined. 


5.2  Predicted  Vahie  Strategy 

When  a  fade  occurs  or  no  taigets  arc  present  in  the  correlation  n^ons,  the  usual 
policy  to  foiiow  is  to  use  the  predicted  coordinate  as  the  measured  one..  This  policy  moves 
the  target  in  a  strai^t  line  trajectory  at  the  last  known  velocity.  This  seems  like  a  sensible 
policy,  since  there  exists  no  a  priori  or  current  information  about  the  status  of  the  target. 

The  effects  of  this  practice  on  the  tracking  equations  are  shown  below.  Equations 
(3.4)  and  (3.5  are  rewritten  as 
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J/fi  =  1000  ft/tM 

i  I  i  I  I  I  !  I  I  lOOOOOOOOOOOOOOOl  I  I  I 
I/f,  =  2000  ft /tec 

IllllOOOOOOOOil 

Pi{.  5.7— 'Selection  and  fade  sequescec  for  fade  zone  Sen^ths  of  40,000  ft  and 
detection  zone  lengtte  of  60.000  ft  fci  two  different  radial  velociiiet; 
r«  4sec 

I/r  =  1000  ft/sec 


I  lOOOOOOOOOOl  10000000000  I  I 

I/p  =  2000  ft/sec 

100000  100000  10 

Fig.  5.8— Detection  and  fade  aeqaencet  for  fade  zone  le.-igths  of  10,000 
ft  and  detection  zone  lengths  of  40.000  ft  for  two  different  radial 
eeiocities;  T  »  4  sec 


.v,j 


.V 


(l--a)  a-c)T 
.  -P/T  (1-/S)  j 

=  [1  r] 


l.V-l 


'J/T\ 


i5.4l 


10.5) 


Equation  f5.5)  is  substituted  into  Eq.  (5.4)  to  obtain 


i%' 


T.V 

- 1 

_ 1 

f — 

o 

{ _ 

1 

i.v-i 


(5.6) 


where 


v'  = 


it  can  be  seen  that  Eq.  (5.6)  is  identical  to  Eq.  (5.4)  but  with  a  and  -cr  equal  to  zero. 
Therefore  Eq.  (5.4)  describes  the  system  undo-  fading  conditions  wh'-»ro  ’ite  parameters  are 
time  varying. 
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5.3  Filter  Noise  Characteristics  Under  Fadii^  Conditions 
The  covariance  Eqs.  (3.18)  and  13.10)  are  rewritten  as 


p’‘'r 

"(1-0)2  2(1  -0)2  (1-0)2 

S-l 

0(2 

1  ! 

\h2\  = 

-;3a-a)  (l-2d>(l-a)  (l-^)(l-a) 

P12 

+ 

i  P02 1 

J 

,32  -2^1-0)  (l-/3)2 

P22 

G.Vp  -  Pli  *  2Pi2  +  P22  •  (5-8) 

The  solutions  for  uq  =  o.Yp/o.Vm  various  fading  sequences,  .say  11001100. are  found 
such  that  if  a  one  occurs,  a  =  0,  p  =  0,  and  if  a  zero  occun,  a  and  take  on  their  ordinary 
values.  The  solutions  are  shewn  in  Figs.  5.9  through  5.13.  Several  general  comments  can 
be  made  concerning  Figs.  5.9  thro'j^  5.13.  First,  tne  heavily  smriothed  filters  (nnall  a 
and  jj)  are  much  less  noisy  than  the  lightly  srucvtl't"!  filters.  In  addition,  the  variances 
seem  to  become  unstable  if  too  many  successive  fades  occur  repeatedly. 

It  was  found  that  the  heavily  damped  filters  were  much  more  susceptible  to  insta¬ 
bilities  than  the  lightly  damped  ones.  A  plausible  explanation  for  this  is  as  follows.  When 
a  fade  occurs  the  estimate  of  the  new  position  can  become  worse  as  time  goes  on.;  WTien 


@  NO  FADES 

(D  f iOOOOOCOOOlOOO--) 
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©  ETC 


unstable  solutions 


Fig.  .S.9— NoTTiiJizwi  predicted  noise  power  js  s  fin.clion  of  tir.x’  f..r 
a  given  fading  '^quence 
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samples 

Fig.  5.1 1— Normalized  predicted  noise  power  as  a  function  Ci  time 
fjT  various  fading  sequences:  fix.  fades,(B)  1000100010  .  <  , 

'©1 1 001 1 001 1 0  ...,(§)  1 1 1 01 1 1 0  .. .  unstable. 

a  =  0.2 

^  °  (A)  NO  FADES 

@(1000100010—) 

©  ( I  1 00 i  1 00 H  0— )  Uf4STABL£  SOLUTIONS 

®  (  mOi  HOI  I  10—) 

®  ETC. 

F.g.  5. 1 2— Normalized  predicted  noise  power  as  a  function  of  time  fc;  various 
fading  sequences:  (^no  fades.C^  1000 1000 10  .< .  .  ,'^11001100110  .  .... , 
'§/  lllOlllOllIO  ..  V  ,(^etc.  The  last  four  are  unstable. 


a  signal  appears  the  heavi'y  damped  system  recovers  very  slowly  and  only  a  poor  estimate 
of  the  velocity  and  new  position  is  found  before  a  lade  again  occurs.  If  too  many  fades 
occur,  tnc  estimate  continues  to  become  worse  The  lightly  smoothed  system  recovers 
more  rapidly  after  a  fade  and  thus  can  handle  lunger  fade  sequences  without  becoming 
unstable.. 

It  can  also  be  observed  that  as  the  number  of  fades  increases  the  noise  power  increases 
rapidly.  One  could  decrease  this  by  sampling  at  a  slower  rate  which  yields  fewer  fades  at 
the  exper.se  of  deterministic  error  increases.  Yet  this  can  lead  to  trouble;  Figs.  5.10  at'id 
5.11  can  be  considered  to  be  sampled  at  two  different  rates.  One  finds  in  Fig.  5.11  that 
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Fig.  5.13~NonnaIized  preditted  noise  power  as  a  function  of  time  for  var¬ 
ious  fading  sequences"  (S)  no  fades,  (S)  110000000000000000001100  .  .  .  , 

©1111 00000000000000001 1 1 100 ... ,®  1 1 11 1 1 000000000000001  mil..., 

©  llllllliOOOOOOOOOOOOllllllll  .  .  .  .©etc.  ®  and  ©  are  unsUble. 


if  the  rate  is  slower  the  system  becomes  unstable  for  shorter  fade  zone  lengths  than  it  would 
have  in  Fig,  5.10.  These  comments  are  not  necessarily  true  when  the  target  moves  through 
several  fade  zones  between  sampling  times.  Alternatively,  rather  than  change  the  sampling 
rate,  could  change  the  target’s  radial  velocity  and  the  same  effect  is  observed. 

If  one  increases  the  sampling  rate,  one  in  general  reduces  a  and  •?.  The  effect  of  in¬ 
creasing  the  lair  pling  rate  with  fades  is  to  increase  the  noise,  and  the  effect  of  reducing 
a  and  /J  is  to  reduce  the  noise.  Unfortunately  the  prediction  noise,  when  fading  occurs,  in¬ 
creases  with  sampling  rate  much  faster  than  can  be  compensated  for  in  changes  in  q  and 
/3  as  can  be  seen  by  comparing  Figs.  5.9  and  5.13.  The  same  effect  would  be  observed  by 
reducing  the  radial  velocity  of  the  target. 


5.4  Deterministic  Errors  Under  Fading  Conditions 

The  error  between  the  predicted  and  the  true  target  position  can  probably  best  be 
studied  by  considering  some  examples.  These  are  shown  in  Figs.  5.14  and  5.15.  Figure 
5.14  shows  that  for  a  realistic  maneuver,  the  deterministic  error  between  the  true  and  the 
predicted  coordinates  was  nearly  6  naut.  mi.  for  about  a  medium-size  fade  zone.  Figure 
r  5.15  shows  that  even  with  fairly  shorl  fade  zones  as  one  might  find  in  the  region  of  30 

naut.  mi.  in  range,  one  can  easily  ho  over  a  mile  in  error  due  to  the  deterministic  process. 
;  If  tl;e  tai^et  turned  such  that  it  flew  parallel  to  the  fade  zone  for  some  time,  the  errors 

could  become  enormou.s.  Of  course  the  total  error  is  larger  than  that  shown  in  these  fig¬ 
ures  due  to  the  addition  of  the  random  errors  discussed  in  Section  5.3. 
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5.5  Probability  of  Breaking  Track 

Recall  from  Section  3.9  that  the  probability  of  breaking  track  depends  on  the  variance 
of  the  predicted  coordinate,  the  deterministic  error  between  Uie  true  and  predicted  coor¬ 
dinates,  and  the  si7e  of  the  correlation  region.  Under  long  fading  conditions  it  was  found 
that  the  variance  of  the  predicted  coordinate  grows  rapidly  and  in  addition  the  determinis¬ 
tic  errors  during  fades  can  become  quite  large.  This  requires  the  correlation  region  to  be 
quite  large  in  order  to  obtain  a  reasonable  probability  of  iireaking  track  for  many  of  the 
pos.sible  fading  sequences  and  target  trajectories.  If  the  correlation  re^on  is  large,  the 
chances  of  tracking  throu^  regions  containing  false  alarms  become  small.  One  could  raise 
the  threshold  to  eliminate  some  of  the  false  alarms  but  then  one  does  not  detect  the  true 
target  as  often,  thus  making  the  fade  problem  worse. 
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5.6  Discussion  of  Results 

The  fading  problem  due  to  the  multipath  effect  was  reviewed  and  a  procedure  for  the 
tracking  equations  to  folio  during  the  time  of  the  fade  was  postulated.  The  variance  of 
the  predicted  posit  '^a  was  found  for  various  fading  sequences  which  depended  for  the  most 
part  on  the  distance  between  the  lobes,  the  radial  velocity,  and  the  sampling  time.  The 
deterministic  error  between  the  predicted  and  true  target  positions  was  found  for  several 
examples.  In  a  general  way,  the  total  error  grows  rapidly  during  a  fade  causing  the  corre¬ 
lation  re^cn  size  to  grow  in  order  to  keep  a  track  established.  Of  course,  other  considera¬ 
tions  such  as  multiple  targets,  clutter,  and  resolution  keep  the  correlation  region  from  be¬ 
coming  too  large. 


6.0  CONCLUSIONS 

The  purpose  of  this  report  was  to  investigate  the  behavior  of  the  a-P  tracker  under 
various  conditions  such  as  maneuvering  targets,  measurement  noise,  false  targets,  and  fad¬ 
ing  conditions.  Well-known  elementary  solution  procedurrc  were  throughout  the 
report. 

'fhe  probability  of  breaking  track  was  found  to  depend  on  the  measurement  noise, 
a,  s-impling  time,  target  trajectory,  correlation  region  size,  false  alarm  structure,  and  fad¬ 
ing  se-quences.  Excluding  fades,  the  system  error  can  be  improved  by  increasing  the  sam¬ 
pling  rate  at  least  until  the  noise  sources  become  correlated.  For  light  filtering  the  noise 
was  higher  but  the  system  tracked  better  in  the  mean.  Conversely  for  heavy  filtering,  the 
noise  was  lower  but  the  error  was  hi^er  in  the  mean.  The  average-target  strategy  seemed 
like  a  good  procedure  to  track  through  light  or  spotty  clutter  as  long  as  the  correlation 
region  was  small  enough  and  the  target  did  not  perform  extreme  maneuvers  during  the 
time  the  false  alarms  were  pre.sent.  The  predicted  value  strategy  for  fades  probably  would 
work  reasonably  well  if  the  target  did  not  maneuver  just  before  or  during  the  ffde,  if  the 
fade  time  were  reasonably  short  com.pared  to  the  detection  zone,  and  if  the  sampling  tim^e 
were  slow  enou^.  If  these  conditions  are  not  met,  quite  large  coiT«lation  legicms  are  re¬ 
quired  which  severely  liir.it  the  ability  to  track  through  clutter  and  to  resolve  multiple 
targets.  Examples  were  given  throughout  the  repo.'t  to  indicate  magnitudes  of  .he  func¬ 
tions  that  are  typically  found. 

Secondary  information  can  also  be  extracted  from  this  report.  The  .-adar  ■,'haracteris- 
tics  can  be  eas  !y  translated  into  tracking  equation  parameters.  For  example,  pilse  widths, 
beam  widths,  scanning  rates,  prfs.  signal-to-noise  ratios,  antenna  pattern  characienstics 
due  to  multipath  effects,  .clutter  rejection  capability,  etc.  can  be  translated  into  measure¬ 
ment  accuracies,  sampling  times,  fade  conditions,  and  false  alarm  characteristics.  Given  a 
particular  radar,  one  should  be  able,  at  lea' :  m  a  rudimentary  way,  to  determine  what  types 
of  tracks  could  be  maintained. 
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